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ABSTRACT 


ie meme ratihy saccepLred that the 10onosphere tilts; 
aay 15 GO Say, an isotonic layer is not at constant 
Poon eaDOvVemmenc SUrtTacecwot the earth. lIonospheric tilt 
has the effect of deflecting a radio ray out of its great- 
circle plane and returning it to earth at an angle not 
Midveol cde wugilemicarine [rom a receiver to a transmitter. 
ine magnitude of error introduced by this effect on radio 
direction finding (RDF) position estimates was studied 
in this paper. A model assigning a tilt bias of less 
Chan three degrees to each RDF station bearing was con- 
structed. Analysis'of a six-station RDF network revealed 
that this amount of tilt has negligible effect on point 


Bmawocmoiml@callon and their confidence regions. 
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I. INTRODUCTION 


Long range radio direction finding (RDF) is the art 
of determining the geographical position of a transmitter 
of unknown location by measuring the azimuthal angle of 
arrival (bearing) of a radio ray at several receiver 
Sites and estimating the location of the transmitter by 
PO rneOLecetaneilacion. RDF techniques typically are 
based on the following two assumptions: 1) The earth is 
modeled sufficiently well by a true sphere, and 2) An 
observed bearing varies from true bearing by random error 
which is distributed normally (gaussian) with mean zero. 
Beeeicdeleis.>o = I + e€, where 0 = observed bearing, 

T = true bearing and e = normal error associated with 
mae receiver. 

Divomuivesuleavulon was Concerned with the possibility 
of the existence of a non-trivial error due to ionospheric 
Peilection Superimposed on true bearing. The model hypo- 
moestizea was 0 = IT + B +t €, where B is a random variable 
which is normally distributed about a non-zero mean. 

This component changes the mean of the distribution of 
PimeouscuvedmOcarinoseiron ! to T + b, where b is the mean 
or 5, 

it 1s well known that reflecting layers in the ionos- 

phere are not of equal height above the surface of the 


mT 2 _* Ad : e e eu wt J | ; A 
Seen ah OT OCces Sn icnosoheric tilt which can provide 





the component B. There are at least two major obstacles 
wommnodcianoe thisee! fect ; 

1) Non-predictable traveling wave disturbances and 
the random phenomena which produce them keep the ionosphere 
idea cOonstany State of flux, and 

2) There does not exist a sufficiently extensive 
Meonmwrorime Meywork te completely map electron densities 
in the ionosphere. These obstacles are discussed in the 
Pequel vogetner with general discussion on the effects 
Of ionospheric tilt on RDF bearings and the estimated 
positions derived from those bearings. A model is pre- 
sired adamconpieuer simulation is wsed toe demonstrate 
Simeon Cumem co ullavcastocatlon and confidence 


moo LoOns . 





Ti. JONOSPHERIC PHENOMENA 


A. GENERAL DISCUSSION 
jne=tenespmaere 2s sa region of electrically charged 

(ionized) air beginning about 25 miles above the surface 
Sieche Garth. ~ Radio Waves Can travel long distances by 
Peace relracuecdnin the TOnospnere and returned to earth. 
It is the refracted ray that is received at the RDF site. 
The degree of ionization and the distribution of the charged 
PTAigeicloom S Mey COnstvane with respect elther to time or 
Beeerapiacal loeavion. iiIihis inconsistency results in 
meommactlOM Ofm@amray in a different manner from one instant 
of time to another for a given Sienal ever a piven patin-. 

| Very grossly the ionosphere may be thought of as a sea; 
macvelo,elayers are not completely flat but contain ripples 
Hike waves on an ocean. In addition to the small-scale 
phenomena that compound inconsistent refraction: 1) First- 
order solar effects and 2) Traveling wave disturbances. 
The nature and effect of these concepts are discussed 


mirst followed by citation of some experimental support. 


B. REMARKS ON PROPAGATION 
Although it is convenient to think of the ionosphere 
wo merrer-like reflecting surface, and such an interpre- 
marlon iS Suiticient for some purposes, in actuality rays 
Mae Dent Or reiracted in the ionosphere before being returned 


PO cartn. The amount of bending or the time a ray spends 


OO 





tHmewneceonmospaere OCefore Deing returned Co earth depends 

in a complicated way on wave frequency, magnetic field 
emdmeaeptn ot penetration. In Figure 1, L represents 
"low-angle" refraction and H "high-angle". It is seen 

that rays arrive at the receiver, R, at correspondingly 

. different elevation (vertical angle). It is also intui- 
mime chat the wnuem—angile ray Spends more time in the ionos-= 
phere since it penetrates deeper. In the SECWiSe ll, IMS eie— 
ence will be made to three layers, E, Fl and F2. Ina 
smooth undisturbed ionosphere, the F layers typically reflect 
Bovgla wien amd ya low ray as in Pigure 1 whereas the E 
Haver typically retleets only one ray. Figure 1 illustrates 
"one-~hop" transmission. Two-hop rays (and higher degree 
hops) result from reflection at the surface of the earth 
back into the ionosphere which again returns the ray to 
PacumagmameGliterent Ceorprapnical location. Terminologi- 
eolly, a Ome—-nop low-angle ray reflected at the Fl region 
Will be designated 1F1L. A two-hop Fe high-angle ray will 
be 2F2H. Similarly, 2E will refer a two-hop ray reflected 
Peewee slayer, For iurther treatment of propagation 
phenomena the reader is referred to any of the number 

Of cextbooks on the subject. Two texts recommended are 
Kelso (1964) and Davies (1965). Also Ames (1964) contains 
Mie xeellene Driet discussion of propagation relating to 


me Crtect of tilt on bearing angle. 
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medium 





Radio ray propagation modes 


FIGURE 1 
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C. FIRST ORDER SOLAR EFFECTS : 

| oan mlntine LoMosphere are leonized by ultra= 
violet rays from the sun and to a less extent by charged 
pareicles Erom the Sul. tnerefore the angle at which the 
Die seo ome Moule tFMestonosphere determines the 
degree of ionization. 

If allowed to ignore the effects of traveling distur- 
Dances, magnetism, earth surface and wind conditions one 
mouldsay thabpeanm tsotomic Layer of the tonosphere is 
highest above the surface of the earth at the equator 
and decreases in altitude with increasing latitude because 
the angle between the sun's rays and local zenith increases 
AolaoeudesiMerecases. | lane result is a north-south Cilt. 
Pmenciieore wT OonuaydirOnm imereases With increasing height 
and higher layers tilt more than lower layers. Figure 2 
exe @eraeca sp Mesperny. Lf the four possible rays reflected 
Polen one eH wlaVers arem=ol mmemameously received at one point 
SrcuVOMlamexpecy, each Fo arrive on a Slightly different 
Teotciieween tore msa sillustraves the difference in eleva- 
tion angles and Figure 3b shows difference in azimuthal 
eagles Of arrival for the four rays. 

PimiddmmioieveombMcmbaukuudanal erfect there exists 
Mecarou-Veoy Uullu. | Blectron density is highest at local 
moon and lowest at local midnight. There is constant 
change throughout the day and changes are very pronounced 
at ionospheric sunrise and sunset. Bramley (1956) estimated 


-~ en” 


the east-west diurnal tilt chenge rate to be on the order 


bd 


dusk 





Great circle plane 


high-angle 
mode 


low-angle | reflected rays 


mode | 


\ \ eS 


Transmitter 
end Receiver 





The effect of tilt on reflection 
of a ray from ionosphere 


FIGURE 2 
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of 0.2 degrees per hour. He does not distinguish between 

Simcoe sana SunSses Wours amd=periods of Bess change. 
Pilot evOomunemGdaltital tilt there exists a seasonal 

Voter ron as the Carth rocates around the sun and the 

sun's direct rays vary between the Tropic of Cancer and 

Mem rOpl cl OmmecmwrrcOrnemmmmnss C€ffect is considered by 

most experimentors to be quite minor and slow to change 

compared to diurnal effects. Munro and Heisler (1963) 

Summarize existing thought on the magnitude and directiona- 


lity of both diurnal and seasonal effects. 


D. TRAVELING WAVE DISTURBANCES 

There are a number of other factors directly influenc- 
mae the shape of an isoionic Paver (Ne toemosonere Gver 
land is considerably different than it is over sea. The 
earth's magnetic field causes drag in the F-layer plasma 
and it varies highly. There are winds, thermal, and coriolis 
Meanee lc One more phenomenon, tCraveling disturbances, 
will be discussed briefly. 

Much literature exists on the subject of traveling 
wave disturbances. Several articles contain compasite 
reviews of existing thought and past experiments and con- 
clusions. One of the best of these is Detert (1965). 

(See also Munro and Heisler (1963) and Heisler (1965)). 

A disturbance is characterized by an increase or 
Gecrease in electron density over background profiles. In 
See xpe iments ouraveling disturbances were included in 
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sophisticated measuring equipment became available they 
became known as traveling wave disturbances and attempts 
HoPeeMacenGeGmmncasure woedir Size and vellocities, and to 
biecote> Chem A wide range of sizes and velocities have 
been reported. Hewish (1951, 1952) reports observing 
-lengths (longitudinal extents) from 2-10 kilometers. 
Bramley (1953, 1955, 1956) observed velocities from 90 
to 1300 km./hr. Chan and Villard (1962) reported dis- 
Minpances Teen L500 km, to over 2000- km. in length and 
velocities from 1450 to 2750 km./hr. Heisler (1963) 
shrewdly points out that size and velocity observations 
depend heavily on the method of measurement. 

Causes of disturbances can be known (e.g., observable 
sun storms) or unknown. Chan and Villard believed the 
large disturbances they observed to have resulted from 
the same event that caused a coincidental change in the 


earth's magnetic field. 


oe GAPE REMENDTAL SUPPORT OF THE TILT CONCEPT 

ine results of three experiments will be presented 
in Support of the existence of a non-zero mean bearing 
error. Sweeney (1970) measured the azimuthal pattern 
realized by a 256-element 2.5 km. broadside array receiv- 
ing HF signals propagated over a 2600 km. east-west path. 
ioawas found that high rays tend to arrive south of low 
rays. soweeney hypothesized that this tendency is a conse- 
meemee Of 19oncespneric tilts havin 
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modeling the ionosphere in a computer ray-tracing program. 
Figure 4 and Figure 5 are reproduced from Sweeney's report. 
The true transmitter bearing was 90 degrees. In Figure 4 
it can be seen that the 1F2L ray oscillates about the 

true bearing during undisturbed conditions. Sweeney sum- 
MerwiZes thay Leow-angle rays are accurate to the resolution 
of his antenna system; that is, one cannot discern a non- 
zero mean component of error in the low-angle rays. Swee- 
ney also concludes that deviations occur mainly in reflec- 
ietietiom tne Garth. if this is in fact the case one 

can draw the following conclusions: 1) deflection out of 
ae cco cle lo anc will nes be measurable in a one—-hop 
low-angle ray, 2) a high-angle ray will be observed noti- 
cably south of the low-angle ray, and 3) multiple-hop 
rays will have more error due to reflections from the earth's 
surface. lit must be kept in mind that Sweeney's two-hop 
rays were being reflected by the Rocky Mountains, an 
Mausuelly rough reflecting surface. 

Bredek (1963) was concerned with round-the-world (RTW) 
propagation. But his comments on the direct ray (Stanford, 
California to Champaign, Illinois) are of interest. 

The direct ray is defined to be the signal received via 

me oeorver Of the two great Circle paths. Bredek observed 
that bearings fluctuated about a daily mean. The winter 
means were north of true bearing and displayed a southerly 


trend until they swung south of true bearing in March. 
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2210 


2220 2230 2240 2250 2300 
Time (GHT) 


Plot of data collected by Sweeney (1970) showing 
higheangle ray arriving south of loweangle ray. 
The data were taken at 23 MHz on 25 August 1969. 
True bearing was 90 degrees, 


FIGURE 4 
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2020 2030 2040 2050 2100 
Time (GMT) 


Plot of data collected by Sweeney (1970) 
illustrating the difference in azimuth of 
different propagation modes, The data were 
taken at 14.5 MHe2 on 30 September 1969. 
True bearing was 90 degrees. 


FIGURE 5 


6 





Whieh crossed true bearing on 21 March, the vernal equinox. 
The conclusion is obvious even if somewhat speculative. 
Mildly stated, bearing means are sensitive to seasonal 
enanee. One otner Bredek observation is worth noting here. 
Variation in RTW bearings was centered "a few degrees" 
i Ot mie sIreCTmbGarines, lhis 1S consistent with 
Sweeney's conclusion. Since the RTW ray has traveled much 
farvhner than the direct ray and is definitely multi-hop 
one expects to observe more variation init. And with 
Sweeney's hypothesis one expects it to arrive south of the 
@irect ray. 

An unpublished experiment conducted at Naval Security 
eeoup ACbivity, okages Island, California, showed the 
existence of an east-west diurnal tilt. Bearings were 
taken at two minute intervals on a signal transmitted from 
Hawaii during ionospheric sunrise and sunset. The data 
Mere analyzea by an autocorrelation function. Over a 
PetwodmOteuWwounours Ghe bearings failed to become statis-— 
tically independent; that is to say, the bearings showed 
meet tm@ve trend to slide in one direction and not fluc-— 
tuate about a cumulative mean. Similar experiments were 
Seneucvced during undisturbed day and night conditions. 

The results were similar to Bramley and Ross (1951). 
Bramley (1953, 1955) and Bain (1955). Bearings showed 
a Siow (up to 20 minutes) quasi-cyclic fluctuation about’ 


a mean in addition to rapid second-to-second fluctuation. 


ine 





eee Ors 1 


There does exist a deterministic element of bearing 
error. At least it can be said that there exists a compo- 
Meme Ol Crromecvhat has a determinable non-zero mean. 
ec Oommen muic com SG bO.GCCeCrmine the effect or this 
non-zero mean component B on RDF "fixes" (estimated trans- 
mitter location) and probability statements about these 
i x CIS : 

According to Dr. Villard” (private communication 30 
November 1971) the present state of technology is such 
ena ane non-zero mean component can be measured to almost 
any accuracy desired. What is lacking is the total commit- 
Rent Tomi ancome Oi DHeMesto Measure SUCH a bias. ~ in 


the absence of this commitment one may account for the 


moron esoy modeling and @statistical methods. 


A. QUALITATIVE FEATURES 

The most severe effect of ionospheric tiit is from 
PoueVesuerile wat 1ON@spheric sunrise and sunset. At other 
times of the day east-west tilt is very gradual. No attempt 
Pili muomccecounuetor sunrise and sunset tilt in this ) 


eden lomsugpesued, however, that this tilt should be 


tor. Cece ard eae. Omeswantord University, is 


Chairman of the Special Committee on Electronics, a panel 
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modeled by an error with positive mean Or bearings an the 
third and fourth quadrants and a negative mean for bearings 
in the first and second quadrants at sunrise and opposite- 
ly at sunset. Between sunrise and sunset (and sunset and 
sunrise) east-west tilt is modeled sufficiently well with 
mean tilt equal to Zero. 

iii mecetoc solemeren—south Cilt 1s more consistent 
than the diurnally periodic east-west tilt. The ionos- 
PHcemeowbilus Wesv Lomeoach in the morning and east to west 
immecvne alvernoon whereas a morth to south tilt exists 
MoCOuectoummemnc ayumi the Sequel error due to tilt, 
peas unau produced by north-south tilt. 

The model is restricted to rays reflected by an ionos- 
phere between sunrise and sunset and to RDF caus and 
targets located in the northern hemisphere. It was con- 
pecveteOPoOmmepresence CONGlttens at mid-latitudes and it 
is assumed to be sufficiently accurate for all latitudes 
ieCcewiciedt is Used ine colleetionser datasfor this thesis. 

The model contains two components of error. One is 
Mc wucliaimiranidem error and the other is due to tilt, B. 
mae tilt component can be reduced to an unbiased random 
component by inserting a zero mean. 

The model assigns a positive mean to error in bearings 
from zero to 180 degrees and a negative mean to bearings 
from 180 to 360 degrees. The error is assumed normal. 
Absolute value of the mean decreases with increasing dis- 
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Uiaeenme Peceiver could not distinguish high-angle rays 
from low-angle rays. As a result the bearing measured 
is due to some unknown combination of high and low-angle 
rays. It was further assumed that the low-angle ray varies 
about the true bearing and the high-angle ray is the sole 
@Cemeributer to error due to tilt. High-angle rays are 
attenuated more rapidly than are low-angle rays. It 
follows that the greater the distance the less effect on 
Dearing will there be due to the high-angle ray and the 
less the deviation due to tilt. 

A systematic standard deviation is assigned to each 
ee stavlon.weluets the basis of the dispersion, of both 
mracoMwcricow cy and tilt error, B. Hor use as a parameter 
in determining e, it increases with increasing distance 
(see Pope (1970)). For use as a parameter of B it decreases 
with increasing distance. Intuitively, the longer a ray 
momen posea LoOrcrror—-producing elements the more dispersion 
ae oom heoecestri bution. This explains the increase 
weecme parameter, call it s, with distance for e. A differ= 
Mireaceumene applies to B. It is claimed that the dis- 
Pao remnoreesis directly proportional to b, the mean of 8. 
The quantity b is in some sense a measure of the strength 
Seeecune eftects that produce error due to tilt. The higher 
mee Valuc taken on by b the more influence on bearing has 
the high-angle ray. Using the same argument as for s 
(i.e., more exposure means higher variability) one con- 


<“— a As pel Fen Ne +L, t- } ) 
Suge s ubas vne higher the b the greater the variance of 
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Eee Datrameter SUpplied in =the model multiplies s by 

b to serve as the standard deviation of B. When b is Zero 
due to great distance (complete attenuation of the high- 
angle ray) s is used as the standard deviation of B. 

A final qualitative feature of the model recognizes 
WicmiaecyS sudabemean error due LO tilt 18 not the same for 
different bearings taken simultaneously from one site and 
mulgems NO Sloped exactly north-south. Predictions of 
actual slope at a given point in the ionosphere are ex- 
Paeeiely Srosseana Quite inappropriate for a given instant 
of time. Also the high-angle ray is mixed with the low- 
Pore rayeiteseme unknown proportion. All of these effects 
mec acceumnued for vin the model by introducing a maximum 
value parameter b' (see Table 1) and multiplying it by a 


Maatorm random variable to produce b. 


B. QUANTITATIVE FEATURES 
The maximum b' was assigned values as a function of 
Geestbance according to Table 1. The values were assigned 
With some uneasiness but an attempt to justify them follows. 
Horedistances less than 50 miles it was assumed that 
a ground wave is predominant and there is no effect from 
ionospheric tilt. At distances greater than 3600 miles 
the high-angle ray was assumed completely dissipated so 
@aat ONCe apain there iS no deflection due to tilt. 
Sweeney (1970) observed high-angle rays that arrived approx- 


itaeely chree degrees south of low-angle rays alone 


JQ 
ey) 
= 


east-west propagation path (see Figure 4 and Figure 5). 
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He modeled this phenomenon with an alpha-Chapman layer in 
the Jones (1966) three-demensional ray-tracing computer 
program using average density data during undisturbed 
eonditions. The results supported the hypothesis that 
the high-angle ray arrives approximately three degrees 
south of the low-angle ray for that particular 90 degree 
- 270 degree path. It is acknowledged that the maxima 
assigned to the distance categories between 900 miles and 
3600 miles are artificial as are the categories themselves. 
Dito cMe asslenmenvts are sufficient to illustrate the 
effect of a component of error due to tilt and they are 


Consistent With the qualitative discussion above. 


24 





HO Oe Unni! Hiro Lo 


ine model was programmed in the Fortran IV language 


Pid cOmpucer Simulation was employed to evaluate the effect 


Sieevioc NOl—ZeroucOmpemenuwol eCrror On Che location of a 
MemOOlnuvaanad wae Sige and Shape Of the confidence regions 
gpenerated by a standard RDF fix technique. Gayobins Ane 

the program were station and target coordinates and a 
Boavlon Systematic Standard deviation. True bearings 
rece COMmpuved and randem and tilt error was Superimposed 
On them. Fix points were computed by a vector method 
(Pope 1971) and the least squares method (Daniels, 1951 
and Kukes and Starik, 1964). Two methods of obtaining 
confidence ere were available, chi-square regions and 
Pivardiace normal regions. Only the latter was used. Two 
random number generators were used to determine bearing 
SietcOr. ss One selected Uniformerandom variates in the inter- 
val (0,1) while the other selected normal random variates 
mer a @iven mean and Standard deviation. Both generators 
pecs uaose recommended by Naylor, Balintfy, Burdick and 
Chu (1966). The basic steps in the simulation are shown 


in Figure 6. 
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Station and target 
coordinates 


Compute distances and 
great-circle bearings 


Add random error 
to bearings 


Add tilt bias 
(mean = 0 for 
"no tilt" runs) 


Calculate initial 
"quick fix" by 
vector method 


Calculate fix by 
least squares method 


Location estimate 


90 percent confidence ellipse 


parameters 





plot 


Basic steps in computer simulation 


FIGURE 6 
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V. SIMULATION 


A network of six RDF stations (Table II) took bearings 
on five targets (Table III). One hundred fixes were com- 
puted for each of the targets assigning only random error 
to the bearings. An additional 100 fixes were computed 
fer Gach target with the model assigning to each bearing 
ECenmooneny Of error due to ionospheric tilt. For each 
set of 100 fixes, means and standard deviations were com- 
puced for latitude and longitude of the f1x point and the 
moanor Semlaxis and minor Ssemiaxis and axis of rotation 
of the 90 percent confidence ellipse. Additionally, the 
fix points were plotted on a graph with rectangular co- 
pmcimaves. —~2his procedure was repeated 25 times for each 


bereceu. 


an 





VI. RESULTS AND ANALYSIS 


The sample means and standard deviations for each si- 
mulation run were tabulated in Tables IV to XIII. Figures 
M@BDCOleMe ic  atce selected Samples Of pious ef fax Locations. 
Meaningful statistical inferences and the validity of the 
model require that the samples collected be random. But 
tie sample mean Latitudes for Omaha, Gimli and Veracruz 
eeow a Cetiawte, trend and overwhelmingly fail the run 
Pio wm rOrtanaommless. Lo 1S interesting to note that the 
mean latitudes for these three targets are the only sets 
of values that fail the test for randomness (with the 
exception of the very stable minor semiaxis values). 

Pren fOr these Targets the mean longitudes show no trend 
whatever. In order to more clearly determine the random- 
ness of the samples of positions, it was observed that 

the data are matched pairs of latitude and longitude. 

The data were reduced to single observations D(i) = lon- 
gitude (i) - latitude (i). The D(i) for the Omaha samples 
pee tcOUlactved in Table XIV. At significance level .05 

the hypothesis that the D(i) constitute a random sample 

is accepted. 

The objective of this thesis has been met without 
further statistical examination. Differences between 
extreme mean locations measure in the low tenths of degrees, 


$ <r on G FI 71 ms & qe ~ ae roa) ean 
Slots eVve se) M1 ices. ine systematic effect of tilt is 


ee 


Micmac att from a practical point of view. 
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One important observation is that the standard devia- 
rons Of the Unmanha axes Of rotation are quite large. 
Although in the Omaha example major and minor semiaxes 
are quite similar in magnitude, a variation from one extreme 
memevlOn amele vO the other involves a displacement of 
epproximately 25 percent of the area of the confidence 


region (crossed area in Figure 13). 


ag 





VII. CONCLUSIONS 


Ka epOlnnamanamc Onl tmeenceureccions Were calculated from 
Pewoeot bearings cConvaining only random error and compared 
tO the points and regions calculated from sets of Dearings 
TomunNmtecm a Component of error due to ilonospherie tilt 
had been added. This component was considered a normal 
random variable. 

ine Ci meeumete superimposing a tilt error on ™bearings 
on the least squares method of computing estimated location 
and confidence regions from a six-station RDF network 
mo Megligible for a méan error due to tilt of less than 
or equal’ Omummece degrees with standard deviation less 
than three degrees. 

The most intriguing development was that only in the 
case where the target was completely surrounded by stations 
did the angular orientation of the confidence region 


womvecapemecoraply., Tilt played no role in this variability. 
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VIII. REMARKS 


The model presented was only a gross approximation to the 
effect of tilt. There are ways available to more accurately 
fame the tonosphere. Predictions of ionospheric character-— 
Poamrcome xis iteeseverale forms. lonespheric Predictions, pub-— 
lished by the Institute for Telecommunications Sciences and 
Aeronomy (ITSA), is a monthly periodical which contains numer- 
teal maps Of Maximum usable frequency at zero range (MUF (ze- 
ro)) and MUF(4000km.) for the F2 layer. In conjunction with 
other publications (e.g., National Bureau of Standards (NBS) 
Circular 462) it is possible to approximate virtual height 
of the F2 layer with respect to geographical location and 
time of day. The NBS Technical Note 40 series is published 
meeeocrly dnameencains predictions of ionospheric electron 
density. These can be used to construct a model based on the 
Chapman layer. (See Barnum (1968) p. 75 and Haydon and Lucas 
(1968)) With these aids the ionosphere can be modeled more 
@ecuravely but in view of the results of this thesis it is 
mipecested tiat an attempt to do so for RDF objectives would 
prove unprofitable. 

The mathematical treatment of the RDF problem is by 
no means new. The theory was well presented by Stansfield 
(1947) and Daniels (1951). Kukes and Starik (1964) present 
amore lengthy and detailed discussion of the same basic 


ee Orye) sary, Kaplan, Keerly, Reeves, and Shaffer 


Syl 





(1966) further defined terms and presented techniques 

for handling the general position finding problem. 
PiotMommeciMinues saseune the earth CO be a true sphere 
oemenploay spherical trlsomemetry, In fact, the earth 

is slightly oblate so that over long distances a correc- 
tion to the spherical treatment is necessary for accurate 
Mecatlen.  In@RBEF @he problem becomes the difference in 
bearing between a spherical great circle and a spheroi- 
dal geodesic. Using parameters for the Clark Spheroid 

of 1866 the difference at mid latitude was found to be 

as great as 0.5 degrees. An additional characteristic 

io uigmee ule Seedesic may Start north of the great circle 
feet Cross CO SOUUmM @S distance increases. In practice 
the detrimental effect of the spherical assumption is 
considered negligible. Thomas (1970) discussed spheriods 
and presented solutions to a geodesic which are adaptable 


PembDoOunamanual and computer calculation. 


ae 





iota Ovo On wun LER SbUDY 


LIGIGUso TMe Cases simulated in this repory were 
PioscimuOnCemoncumave UNE "eirect of Gilt on RDF Fix points, 
Mae miiveoresuinewslde—ci feCe provided an outrall., Whe 
Salentatdon of a confidence ellipse can change consider-— 
piely Coe ceorrapnicale areca covered by that région. Fur— 
icunioure lL G mame catenuUmMave one Variability of orlentation 
mewocras eM the Lecation of the tarcet relative to the 
RDF network. The variability of orientation of confidence 
merions Can have tremendous impact on the validity of 
probability statements based on RDF techniques. A study 
of the effect on confidence regions of target location 
and network configuration is suggested. The orientation 
minonwcommacenee recion depends on the magnitude of Wariance 
Meecachuwcarimno used to calculate the fix point. A study 
momaebermmune tie Sensitivity of this orientation as a 
mance lon Of variance in individual bearings may prove 


valuable. 
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APPENDIX A 


TABLE 1 


Distance 
(miles) 


0~50 
50-900 
900-1800 


1800-2700 


2700-3600 


> 3600 


b? 
Maxima 

(degrees) 
0.0 
320 
2.5 
2.0 
1.5 
0.0 





Maxima assigned for the calculation of 
means for the component of bearing error 
due to ionospheric tilt as a function 


of distance. 
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TABLE II 


Station Location 

Lat. Lon. 
Miami 25.6 80.2 
Boston 42.4 71.0 
Grand Falls 48.0 94.0 
Seattle 47.5 122.5 
Los Angeles ~ 34.0 118.5 
Houston 30.0 97.9 


Composition of RDF network for simulation 


TABIE III 


Location 
at. Lon. 


Omaha 
Veracruz 
Gimli 
Atlantic 


Pacific 





a 


Tarrets 8EmN * : a) 1} {4 ry 
Ge tiv vO Wwe VA in a? ated U a Or 





TABLE IV 


Location 90 Percent Confidence Ellinvse 
rinor Axis of 


Semiaxis SD Rotation 


a jor 


1.295 
1.296 
1.295 
1.296 
1.294 
1 $294 
1.295 
1.295 
1.296 
1.296 
e2o7 
1.294 
1.294 
1.295 
1.296 
1.295 
12295 
1.295 
1.297 
1.295 
1.296 
1.295 
1.296 
1.293 
1.295 


SDs |Semilaxis. SD 


Sample means and standard deviations 


Targot: 


Omaha 


@ e ¢ 
oOoooqoqooqo°o 


BIGGRRAS 


oOo99090 [2909 
e ¢ 


113124 
10.742 
10.454 
10.343 
10.948 
10.517 
10.253 
10.562 
10.811 
11.686 
11.108 
11.106 
10.027 
11.190 
10.82 
10.445 

9.985 
11.283 

9.916 
10.046 
11.090 
11.408 

9.772 
10.921 





Mode: No tilt 


All values in degrees 


seh 





TABLE V 


Location 90 Percent Confidence Ellinvss 
Major Minor Axis of 
SD |Isemieaxis SD Semiaxis Sp Rotation 





41.112 


Sample means and standard deviations 
Target: Omaha 
Mode : TLC 


All values in degrees 
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TABLE VI 


Location 90 Percent Confidence Ellivse 
Major Minor Axis of 
SD |(Semiaxis SD Semiaxis SD Rotation 





572968 


Sample means and standard deviations 
Target: Atlantic 
Mode: No tilt 


All values in degrees 
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TABLE VII 





Location 90 Percent Confidence Ellinse 
Major Minor Axis of 
SD Lon. SD jSemiexis SD Semiaxis Sb Rotation SD 
0.085 |!0.901 O.001 | =26-142 |Oe210 
0.072 {0.901 0.001 | ~26.030 | 0.211 
0.071 |0.901 0.001 | ~26.070 | 0.197 
0.086 10.901 | 0.001 | -26.103 | 0.209 
0.090 |0.901 | 0.001 | -26.168 | 0.228 
0.085 |0.901 0.001 | =26.149 | 0.232 
0.068 [0.901 0.001 | =26.220 | 0.254 
0.075 {0.901 0.001 | +26.137 | 0.220 
0.096 {0.901 0.001 | -26.220 | 0.211 
0.088 [0.901 |0.001 | -26.178 | 0.243 
0.078 10.901 0.001 | #26.165 | 0.246 
0.067 (0.901 {0.001 | -26.175 10.199 
0.074% 10.901 0.001 | -26.241 10.205 
0.071 10.901. |0.001 | ~26.200 |0.207 
0.079 0.901 0.001 | -26.201 |0.221 
0.068 0.901 0.001 | -26.107 |0.184 
0.098 0.901 {0.001 | -26.128 {0.251 
0.082 0.901 0.001 | -26.139 |0.180 
0.084 0.901 0,001 | 26.131 10.215 
0.083 0.901 0.001 | -26,.088 |0.223 
D.088 0.901 0.001 | -26.186 |0.224 
0-076 D.901 0.001 | -26.183 10.168 
D.081 P.901 0.001 j -26.163 10.201 
D.072 0.901 0.001 | -26.192 |0.223 
0.083 p.901 0.001 | -26.214 |0.225 


Sample moans and standard deviations 


Target: 
Mode: 


Atlantic 
Tale 


All values in degreos 
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Location 


SD 


0.152 
0.154 
0.181 
0.163 
0.183 
0.179 
0.148 
0.169 
0.160 
0.160 
0.143 
0.190 
0.179 
0.168 
0.186 
0.170 
0.168 
0.179 
0.203 
0.185 
0.203 
0.193 
0.194 
0.157 
0.144 


Lon. 


135.044 
135.140 
134.981 
134.889 
135.083 
135.034 
135.001 
135.096 
134.884 
135.008 
135.016 
134.919 
135.110 
134.977 
134.906 
135.006 
135.028 
135.054 
135.030 
135.053 
135.090 
135.035 
134.998 
135.017 
135.032 


TABLE VIII 


oD 


0.502 
0.580 
0.703 
0.585 
0.641 
0.631 
0.522 
0.614 
0.562 
0.557 
0.549 
0.674 
0.639 
0.645 
0.647 
10.587 
0.625 
0.642 
0.754 
0.670 
Oey2t 
0.677 
0.616 
0.631 
0.583 





Major 


Semiaxis 


Minor 


SD Semiaxis 


Samplo means and standard deviations 


90 Percent Confidence Ellinse 


Axis of 


SDaeekeGvation 


20.651 
20.721 
20.747 
20.658 
20.639 
20.633 
20.788 
20.692 
20.697 
20.699 
Z0e7o9 
20.723 
20.748 
20.674 
20.796 
20.654 
20.697 
20.693 
205755 
20.684 
20.689 
20.702 
20.701 
20.739 
20.652 


SD 


Target: Pacific 
Mode: No tilt 
The dimension of all values is degrees 
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TABLE IX 


Location 90 Percent Confidence Ellinse 
Axis of 
SD Rotation 


Major Minor 
Semiaxis SD Semiaxis 


0.050 | 0.906 
2064 | 0.906 
2054 | 0.906 
~072 | 0.906 
-073 | 0.906 
0.067 | 0.906 
0.080 | 0.906 
0.069 | 0.906 
0.067 | 0.906 
L077 0.906. 
0.070 | 0.906 
0.076 | 0.906 
0.069 | 0.906 
0.082 | 0.906 
D.069 | 0.906 
0.073 | 0.906 
~072 |0.905 
0.057 | 0.906 
0.080 10.906 
0.072 19.906 
D064 10.906 
0.080 | 0.906 
D074 | 0.906. 
HOOT wel tO coco 
0.077 | 0.906 





Sample means and standard deviations 
Target: Pacific 
Mode : Tilt 


All values in degrees 
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oD 











Location 


Sb Lon. 


96.970 
$6 ..975 
96.977 
6.979 
$6.962 
$5 .977 
6.970 
96.953 
95.977 
0% S58 
$6 ..975 
05.981 
G5 988 
95 .967 
S56 .562 
96.954 
96.952 
GO 952 
65.975 
96.957 
SO .S49 
96.930 
65.986 
05.950 


96.985 


TABLE X 


90 Percent Confidence Ellinse 


Major 


SD iSemisxis 


0.053 
0.072 
0.051 
C253 
Oe 


1.415 
1.415 
1.418 
1.418 
1.417 
1.417 
1.418 
1.419 
1.419 
1.418 
1.420 
1.420 
1.421 
1.422 
1.421 
1.416 
1.416 
1.419 
1.418 
1.419 
1.419 
1.419 
1.418 
1.420 
1.416 


SD Seniaxis 


0.008 
0.010 
0.010 
0.009 
0.009 
0.008 
Oe ee 
0.010 
0.012 
0.010 
orien ae) 
0.010 
Oe 007 
Olan 
0.011 
0.008 
0.008 
0.011 
0.009 
0.011 
0.010 
0.009 
0.003 
0.010 
0.009 


flinor 


1.416 
ake 
es ke 
1.114 
1,115 
eis 
ee ey 
dopach 
Teese 
oka 
PRs 
pe Eolas 
1.133 
1, dol 
pte ales: 
pea 
Lea 5 
eseee 
eee 
eee 
oie 
1.114 
1.114 
Les 
1.115 


Sample means and standard ceviations 


Targets 
Mode! 


Gimnly 


No tilt 


All values in dcogrees 
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0.004 
0.005 
0.005 
0.004 
0.005 
0.004 
0.006 
0.005 
0.006 
0.005 
0.005 
0.005 
0.004 
0.905 
0.005 
0.004 
0,004 
0.005 
0.0C4 
0.005 
0.005 
0.004 


0.0C4 


0.005 
0.004 





Axis of 


SU Kotatiaon 
234967 
34,088 
B35500 


SSPV ae) 


23.050 
33.703 


33.489 


Bones 


Boe sor 
895257, 
6 sie olen 
33.439 
Be s7ee 
32 6943 
33.149 
Soro 
33.E42 
Sees 
33.050 
Bins 
BBeece 
33.467 
33.447 


33309 


go.205 


SD 















































































TABLE XI 


90 Percent Confidence Ellinss 





Target: 


Mode : 


ye. 


Giml4 
Tiit 


All values in degrees 


Sample means and standard deviations 


Ita jor Minor Axis of 
Spemawts. ove cOniaxis. SU “Rotation SD 
1.420 | 0.010] 1.113 | 0.005] 30.952 11.273 
fie eeoll | tlie 10.005 391.338 idee 
i220) 102007 | 1.81 3 1-0, 003) all.276 ioe 
1.420 | 0.009 | 1.113 |0.0C4 | 31.443 11.076 
Perle OROC7 | 1.itt 10,00) sites  simmesg 
.1.419 | 0.005 | 1.114 | 0.005 | 31.185 {1.050 
S20 eOm0C9 | 1e113010,0C4 | 30ce73 Biamee2 
eG GeOLO | Walls 10.005 i130.720. alec er 
e220 GOnOCG 4 e113 4 0.0CH i 20.712) a bs002 
4.420} 0.0101 1.114 10.005] 30.696 {1.154 
Pee lOmmomole lel. irr O, OCS | 30.510) wei Ge 
les Peden tool 2 leet 3.170 .005)| sOn5so meek S 
1.419 | 0.009 } 1.411% | 0.00% | 30.440 [1.038 
1.419 | 0.009 | 1.114 10.005 | 30.629 {0.754 
I N2OROeOLO len t1 3 0,005 | 262179 mie. Cos 
Peo OmperOLO iml 10.005 | Gine7 3 mies 
iS 20eO sol Ouest ome. 005 |) 3076 melenc-e 
1.418 | 0.010 | 1.114 | 0.005 | 31.198 {1.133 
eT OMmOsOLO fides ite 10.005 | 315175 meen ol 
ie eon Ol Ollie 14 90.005 | 30.G0Gumiia oct 
iS Os009 Wleel 149 0.00% | 315053 mltni70 
1.419 | 0.008 |1.114 |0.0C4 | 30.950 {0.910 
te USO OOS edi) | 0,004 | 30. 91e amleeie7 
eeu Os Ol 2 eet SMD. OCS | 30.767 milenere 
14522740011 Weeeat2iG.005 | 30.488 941.065 





TABLE XII 


Location 90 Percent Confidence Ellivse 
Major Minor Axis of 
Semiexis SD Semiaxis SD Rotation 





jie 
e e e e e e e e 


RARBRARARSSSASAASSSAAAAABR 


e 
Wau 


RAAK 


CON OMONFRPNNW OW 


e ® ® t @ o ° e 
e © e e e e @ ® e 


SEGESELLESES 


oe 
® 


€ 
WN 
WwW 


Un in ON ON a nn 


® @ 


fab eb Fe PD pe Fe pb pe pb pd FD eb fd FD pd pd pd pe fd pd pd Fd Fd PD fs 


1 
1 
1 
1 
1 
1 
1 
1 
i 
1 
i 
i 
1 
1 
iL 
i 
1 
1 
2 


PERE 


e 


95-996 


Sample means and standard deviations 
Target: Veracruz 
Mode : No tilt 


All values in degrees 


at 








18.472 


Location 





Lon. 


95.258 
95.872 
95.858 
95.934 
95.952 
95.882 
95.963 
95.44 
95-909 
95.881 
95.904 
95.874 
95.887 
95.918 
95.927 
95.910 
95.970 
95.245 
95.935 
95.045 
95.899 
95.926 
95.518 
95-905 
95.968 


ee ES LT TT Tae * 


TABLE XIII 


oD 


Ueceu 
0.245 
On 258 
0.150 
0.126 
0.254 
0.073 
0.270 
0.195 
0.272 
0.375 
0.271 
0.239 
0.284 
0.169 
0.191 
0.052 
0.292 
0.152 
0.205 
0.222 
0.154 
One? 
0.204 
0.095 


ba jos 


= * 
Semicexis 


oOo 
Ae? 


O Percent Confidence Ellinse 


SO soem iay1s 


0.011 
0.009 
0.009 
0.010 
0.009 
0.009 
0.910 
0.011 
0.009 
O20l2 
0.010 
g2010 
0.003 
0.009 
O7erl 


Minor 


1.057 
1.058 
1.058 
1.058 
1.058 
1.058 
1.057 
1.053 
1057 
De 57 
1.053 
1.057 
1.057 
15057 
1.057 
1.059 
1.059 
1.058 
1.058 
1.058 
1.058 
SOY, 
1.053 
1.058 
1.057 


Sample means and standard deviations 


0.003 


Axis of 
SU heprataon 


~£6 989 
~87 .139 
86 917 
~87 .028 
~-57 .097 
~&7 .129 
~87 .240 
@57 037 
=-27.180 
#87 .211 
-87.073 
-37.189 
-87 361 
-57 2599 
=C/7 77 
~85 6989 
-87.019 
~87 135 
~87 059 
~87 157 
~87.195 
=-87.072 
~87 C24 
-87 2089 
=87 «240 











Veracruz 
ible, 


Targot: 
Mode: 


All values in degrees 


Die 





TABLE XIV 





Reduced location data 
for the Omaha samples 


5 
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